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Abstract 


Results of a search for new phenomena in events with large missing transverse momentum and a 
Higgs boson decaying to two photons are reported. Data from proton-proton collisions at a center-of- 
mass energy of 8 TeV and corresponding to an integrated luminosity of 20.3 fb _1 have been collected 
with the ATLAS detector at the LHC. The observed data are well described by the expected Standard 
Model backgrounds. Upper limits on the cross section of events with large missing transverse mo¬ 
mentum and a Higgs boson candidate are also placed. Exclusion limits are presented for models of 
physics beyond the Standard Model featuring dark-matter candidates. 
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Although the existence of dark matter (DM) is well 
established, nearly nothing is known of its underlying 
particle nature [T] . Many DM candidates have been pro¬ 
posed, and attempts made to connect them to physics 
beyond the Standard Model (SM) at the scale of elec- 
troweak symmetry breaking [5] that would naturally ac¬ 
commodate the observed relic density [3]. 

Collider searches for weakly interacting dark matter 
rely on the inferred observation of missing transverse mo¬ 
mentum 0J E™ lss recoiling against a visible final-state 
object X, which may be a hadronic jet 0 0, photon 
( 7 ) [21 [ 8 ], or W/Z boson [MI]. The discovery of a Higgs 
boson [I21H3] (H) creates a new opportunity to search for 
beyond-the-SM (BSM) physics giving rise to H + E™ lss 
signatures on ng. In contrast to the aforementioned 
probes, the visible H boson is unlikely to be radiated 
from an initial-state quark or gluon. This has the impor¬ 
tant consequence that the H + E™ lss signature directly 
probes the structure of the effective DM-SM coupling; 
see Fig. [Tj 

If the mass of the DM particle is less than half of the 
Higgs boson mass ttih, the Higgs boson may decay di¬ 
rectly to DM. Such decays have been searched for us¬ 
ing LHC data, and null results provide powerful con¬ 
straints on the invisible branching ratio of the Higgs bo¬ 
son in several different production modes including WH 
or ZH dunlin], and qqH da ng. However, the mass 
of the DM particle may be larger than mjy/2, in which 
case these searches are not sensitive, and approaches such 
as analysis of H + E™ lss events are required. 

Two approaches are commonly used to model generic 
processes yielding a final state with a particle X recoiling 
against a system of noninteracting particles. One option 
is to use nonrenormalizable operators in an effective field 
theory (EFT), which is agnostic about the details of the 
theory at energies beyond the experimental sensitivity. 
Alternatively, simplified models that explicitly include 
the particles at higher masses can be used. The EFT ap- 



FIG. 1: Schematic diagram for production of DM particles \ 
in association with a Higgs boson in pp collisions, mediated by 
electroweak bosons ( H , Z , 7 ) or new mediator particles such 
as a Z' or scalar singlet S. The gray circle denotes an effective 
interaction between DM, the Higgs boson, and other states. 


proach is more model-independent, but is not valid when 
the typical momentum transfer approaches the scale of 
the high-mass particles that have been integrated out. 
Simplified models do not suffer from these concerns, but 
include more assumptions by design and are therefore 
less generic. The two approaches are thus complemen¬ 
tary and both are considered here. 

In this Letter, results are reported from a search for 
H + E™ lss events in data collected by the ATLAS detec¬ 
tor from pp collisions with center-of-mass energy y/s = 
8 TeV and corresponding to an integrated luminosity 
of 20.3 fb~ 1 , produced by the Large Hadron Collider. 
The H 77 decay mode is used exclusively, as the 
small branching ratio is mitigated by the distinct dipho¬ 
ton resonance signature and the low expected number 
of background events with significant E™' ss |14) . AT¬ 
LAS measured previously the differential cross section of 
H —)► 77 production with respect to several kinematic 
quantities [ 20 ], including E™ lss ; the search reported here 
uses a subset of those data optimized for sensitivity to 
production of dark matter in association with a Higgs 
boson. 

The ATLAS detector (21. is a multipurpose particle 
physics experiment with a forward-backward symmetric 
cylindrical geometry and nearly 47 r coverage in solid an- 
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gle. Events were selected using a trigger that requires two 
photons, with leading (subleading) E T > 35 (25) GeV. 

A photon is reconstructed as a cluster of energy with 
\r/\ < 2.37 deposited in the electromagnetic calorimeter, 
excluding the poorly instrumented region rj £ [1.37,1.56]. 
Clusters without matching tracks are classified as uncon¬ 
verted photon candidates. The photon energy is cor¬ 
rected by applying an energy calibration derived from 
Z — > e + e~ decays in data and cross-checked with J/^ —► 
e + e“ and Z —► decays in data [21] ■ Identification 

requirements are applied in order to reduce the contam¬ 
ination dominantly from 7r° or other neutral hadrons 
decaying to two photons. The photon identification is 
based on the profile of the energy deposit in the first and 
second layers of the electromagnetic calorimeter. Pho¬ 
tons have to satisfy the ‘tight’ identification criteria of 
Ref. [23] . They are also required to be isolated, i.e. 
the energy in the calorimeters in a cone of size A R = 
sj (A77) 2 + (A (j>) 2 = 0.4 around the cluster barycenter, 
excluding the energy associated with the photon cluster, 
is required to be less than 6 GeV. This in-cone energy 
is corrected for the leakage of the photon energy and 
for the effects of multiple pp interactions in the same or 
neighboring bunch crossings superimposed on the hard 
physics process (referred to as pileup interactions) [24] , 
Finally, for each photon the scalar sum of the transverse 
momenta pt of tracks originating from the diphoton ver¬ 
tex with pt > 1 GeV and Al?(track,cluster) < 0.2 must 
be less than 2.6 GeV. The diphoton production vertex is 
selected from the reconstructed collision vertices using a 
neural-network algorithm as described in Ref. [ 23] . 

The momentum imbalance in the transverse plane is 
obtained from the negative vector sum of the recon¬ 
structed and calibrated electrons, muons, photons and 
jets and is referred to as missing transverse momen¬ 
tum _E“ 1SS . The symbol E™ lss is used for its magni¬ 
tude. Calorimeter energy deposits are associated with 
a reconstructed and identified high-px object in a spe¬ 
cific order: photons with px > 10 GeV, electrons with 
Px > 10 GeV, and jets with p T > 20 GeV. Deposits not 
associated with any such objects are also taken into ac¬ 
count in the -Ex‘ ss calculation [25] using an energy-flow 
algorithm that considers calorimeter energy deposits as 
well as inner-detector tracks [26] . The energy resolution 
is typically 11% near the threshold at 100 GeV for the 
considered signal scenarios. 

Quality requirements are applied to photon candidates 
in order to reject those arising from instrumental prob¬ 
lems. In addition, quality requirements are applied in 
order to remove jets arising from detector noise or out- 
of-time energy deposits in the calorimeter from cosmic 
rays or other noncollision processes m- 

Selected events are required to have a Higgs boson can¬ 
didate consisting of two photons with diphoton invariant 
mass m 77 € [105,160] GeV with transverse momenta 
satisfying leading (subleading) p 7 > 0.35(0.25)to 77 . In 


addition, large missing transverse momentum is required, 
E™ lss > 90 GeV, as well as large transverse momentum of 
the 77 system, pP > 90 GeV in order to suppress back¬ 
ground events where E™ lss is caused by mismeasurement 
of the energies of identified physics objects. These se¬ 
lection requirements were derived by optimizing the ex¬ 
pected upper limits on H + Ep ss production for the set 
of models described below. 

Contributions to the 77 + F™ lss sample from SM pro¬ 
cesses include those that produce a Higgs boson in as¬ 
sociation with undetected particles (predominantly ZH 
with Z —^ vv and WH with W £v) as well as non¬ 
resonant diphoton production ( 77 , W 77 , -Z 77 ), W 7 and 
Z^j production where an electron is misidentified as a 
photon, and photon+jet production in which the jet is 
misidentified as a photon. 

Samples of simulated events are used in order to 
measure the efficiency of the selection for dark-matter 
models, as well as to estimate the contribution of SM 
H + E™ lss processes. Contributions from other back¬ 
ground processes are estimated from ?n 77 sidebands in 
the data. 

Following the notation of Ref. [14] , a set of EFT mod¬ 
els are considered in which the effective operator La- 
grangian term can be written as |x| 2 |i7| 2 , X^l 5 X\H\ 2 , 
X*d^xH 1 D/jH, or xi^xB^H^ D u H : where the DM field 
X is a scalar in the first case and a fermion in the remain¬ 
ing cases and B ^ is the U(1)y field strength tensor. The 
interactions of SM and DAI particles are described by two 
parameters: the DM particle mass m x and the suppres¬ 
sion scale A of the heavy mediator that is integrated out 
of the EFT. In a theory that is valid to arbitrary energies 
(ultraviolet complete), the contact interaction would be 
replaced by an interaction via an explicit mediator V. 

In addition, simplified models El with a massive 
vector ( Z '), or a scalar ( S ) intermediate boson are 
tested. All H + Ep ss DM models are generated 
with Madgraph5 [2B] version 1.4. 8 .4, with shower¬ 
ing and hadronization modeled with Pythia8 [22: ver¬ 
sion 1.6.5 using the AU2 parameter settings E; the 
MSTW2008LO [31] parton distribution function (PDF) 
set is used. Values of m x from 1 to 1000 GeV are con¬ 
sidered. Production of ZH and WH is modeled with 
Pythia8 using CTEQ6L1 PDFs [52] . Samples are nor¬ 
malized to cross sections for WH and ZH production cal¬ 
culated at next-to-leading order (NLO) [33], and next-to- 
next-to-leading order (NNLO) [Ml hr QCD, respectively, 
with NLO electroweak corrections (35] in both cases. 

Differing pileup conditions as a function of the instan¬ 
taneous luminosity are taken into account by overlaying 
simulated minimum-bias events generated with Pythia8 
onto the hard-scattering process such that the observed 
distribution of the average number of interactions per 
bunch crossing is reproduced. The simulated samples 
are processed with a full ATLAS detector simulation [36] 
based on Geant4 m and a simulation of the trigger 
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FIG. 2: The best-fit background estimates to the 18 observed 
events are 14.2 ± 4.0 (continuum backgrounds) 1.1 ± 0.1 (SM 
Higgs boson backgrounds) and 2.7 ± 2.2 (BSM Higgs bo¬ 
son), including both statistical and systematic uncertainties. 
An unbinned maximum-likelihood fit to the spectrum is used 
to estimate the number of events from the continuum back¬ 
ground and from H —> 77 decays; the individual components 
are shown as well as their sum. 


system. 

To distinguish contributions from processes that in¬ 
clude H —> 77 decays from those that contribute to the 
continuum background, a localized excess of events is 
searched for in the to 77 spectrum near the Higgs boson 
mass, run = 125.4 GeV. Probability distribution func¬ 
tions that describe the H 77 resonance or the contin¬ 
uum background are defined in the range 105-160 GeV as 
described below. The contributions from each source are 
then estimated using an unbinned maximum-likelihood 
fit to the observed m 77 spectrum. 

The m 77 spectra of the signal models of .ff+DM pro¬ 
duction and SM Higgs boson background processes are 
modeled with a double-sided Crystal Ball [35] function; 
the width and peak positions are fixed to values extracted 
from fits to simulated samples. An exponential function, 
e am 11 w ith free parameter a is used to describe the to 77 
distribution of the continuum background. The chosen 
continuum fit function is validated using simulated sam¬ 
ples of the irreducible background processes and in three 
data samples adjacent to the signal region, but with re¬ 
laxed requirements on U™ lss , on p^ 7 , or on photon iden¬ 
tification. Results of the fit to data in the signal region 
are shown in Fig. [2j 

Systematic uncertainties from various sources affect 
the number of SM Higgs boson events in the resonant 
background, the predicted shape and location of its peak, 
as well as the efficiency of the selection for the signal 
models considered. 

The uncertainty on the integrated luminosity, 2.8%, is 
derived following the same methodology as that detailed 
in Ref. [39] using beam-separation scans. Uncertainties 
on the efficiency of the photon isolation requirement, pho¬ 


ton identification requirement, and trigger selection are 
measured in an inclusive SM Higgs boson sample to be 
2.8%, 2.1%, and 0.2%, respectively. Uncertainties in the 
photon energy scale and resolution lead to respective un¬ 
certainties of 11% and 0.3% in the position and width of 
the H 77 peak. Additional uncertainties on the jet 
energy scale and resolution as well as the calibration of 
unclustered hadronic recoil energy contribute to uncer¬ 
tainty in the E™ lss , leading to 1 . 2 % uncertainty on the 
efficiency of the selection for the signal models from the 
U™ lss and requirements. The impacts on the selec¬ 
tion efficiency of the uncertainties on the levels of initial- 
state and final-state radiation are assessed by varying the 
Pythia8 parameters, as in Ref. these are found to 
be typically at the level of 1%. The total uncertainty on 
the selection efficiency for peaking SM Higgs backgrounds 
and signal models is 4.0%. 

The theoretical uncertainties on the WH and ZH pro¬ 
duction cross sections come from varying the renormal¬ 
ization and factorization scales and from uncertainties 
on the parton distribution functions eh HDHia follow¬ 
ing the PDF4LHC prescription. The Higgs boson decay 
branching fractions are taken from Refs. 431144] and their 
uncertainties from Refs. [JSJ 0B]. The total theoretical 
uncertainty on the H + E™ lss contribution is 6 %. 

The number of events observed in the data corresponds 
to a 1.4 a deviation using the asymptotic formulae in 
Ref. (47]. As the events observed do not include a sta¬ 
tistically significant BSM component, the results are in¬ 
terpreted in terms of exclusions on models that would 
produce an excess of El + U™ ISS events. Upper bounds, 
detailed below, are calculated using a one-sided profile 
likelihood ratio and the CL$ technique [48] 09], evalu¬ 
ated using the asymptotic approximation m, which was 
ensured to be valid for the available number of events. 

The most model-independent limits are those on the 
fiducial cross section of H + E™ lss events, including SM 
and BSM components, ax A, where cr is the cross section 
and A is the fiducial acceptance. The latter is defined 
using a selection identical to that defining the signal re¬ 
gion but applied at particle level, where £!“ lss is the vec¬ 
tor sum of the momenta of the noninteracting particles, 
photon isolation requirements are not applied, and a sim¬ 
pler requirement on photon pseudorapidity \rj\ < 2.37 is 
made. The limit on cr x A is derived from a limit on the 
visible cross section ax Axe, where e is the reconstruction 
efficiency in the fiducial region. An estimate e = 56% is 
computed using the simulated signal samples described 
above with no quark or gluon produced from the main 
interaction vertex; the efficiencies vary across the set of 
models by less than 10%. The observed (expected) up¬ 
per limit on the fiducial cross section is 0.70 (0.43) fb at 
95% confidence level (CL). These limits are applicable to 
any model that predicts H + E!f lss events in the fiducial 
region and has similar reconstruction efficiency e. 

Limits on specific models of BSM H + E™ lss produc- 
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FIG. 3: Profile likelihood ratio (A) as a function of OBSM.fid, 
the fiducial cross section for production of a BSM H+DM 
process in the 77 + E™ lss channel taking into account the 
contribution of the SM component. The solid blue likelihood 
curve shows that the number of events observed in the data 
corresponds to a 1.4 tr deviation using the asymptotic for¬ 
mulae in Ref. m- The dotted green likelihood curve only 
includes statistical uncertainties. The dashed red likelihood 
curve allows for modifications of the central value and uncer¬ 
tainty on the SM component as described in the text. 

tion depend on the prediction of the H+E™ lss component 
produced via ZH or WH; calculations of this theoreti¬ 
cal quantity will improve with time and may depend on 
the details of a specific BSM theory. Following the pro¬ 
posal of Ref. [50], the profile likelihood ratio of the cross 
section for BSM 7L+DM production in the 77 + E™ lss 
channel is provided with the SM component fixed to 
the central value of the theoretical calculation, which al¬ 
lows later reinterpretation for any modified prediction 
and uncertainty, as shown in Fig. [3] This approach re¬ 
quires knowing how a change in the SM-like component 
modifies the best-fit BSM component; in this case where 
the SM-like and BSM components are indistinguishable, 
A-^bsm = — ATVgM-iike- The limits on the parameters 
of the specific BSM models considered in this Letter are 
calculated using the prediction and uncertainty for the 
SM component as described above. 

Limits on DM production are derived from the cross- 
section limits at a given DM mass m x , and expressed 
as 95% CL limits on the suppression scale A or coupling 
parameter A for the effective field theory operators; see 
Fig.[4]for limits for D^H and xi^xB^H^D v H 

operators. For the lowest m x region not excluded by 
results from searches for invisible Higgs boson decays 
near m x = mu/ 2, values of A up to 6 , 60, and 150 
GeV are excluded for the X* 75 Xl-H| 2 , X ( d^X^'D^H. 
and xi^xBD v H operators, respectively; values of 
A above 25.6 are excluded for the |x| 2 |-H | 2 operator. As 
discussed above, the effective field theory model becomes 
a poor approximation of an ultraviolet-complete model 




FIG. 4: Limits at 95% CL on the mass scale A as a function of 
the DM mass ( m x ) for two of the four EFT models considered. 
Solid black lines are due to H + E™ las (this Letter); results 
where EFT truncation is applied are also shown, assuming 
coupling values g = ^Jg q g x = 1, 47 r. The g = An case overlaps 
with the no-truncation result. The blue line indicates regions 
that fail the perturbativity requirement of g < Air, the red 
line indicates regions excluded by Z boson limits m on the 
invisible branching fraction (BF), and the pink line indicates 
regions excluded by the LUX Collaboration EH- 


containing a heavy mediator V when the momentum 
transferred in the interaction, Q tr , is comparable to the 
mass of the intermediate state my = A y/g q g x |54j, ^5], 
where g q and g x represent the coupling of V to SM and 
DM particles, respectively. To give an indication of the 
impact of the unknown ultraviolet details of the theory, 
limits are computed in which only simulated events with 
Qtr = 'm'xx % m v are retained; these limits are shown 
for values of x /g q g x = 1 or 47 r in Fig. |4| This proce¬ 
dure is referred to as truncation. In addition, limits are 
derived on coupling parameters for simplified models as 
shown in Fig. [5] For a vector-mediated model, limits 
are placed on the coupling g q of the mediator to quarks, 
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FIG. 5: Limits on coupling parameters for simplified models 
with a heavy mediator with mass of 1 TeV. All constraint 
contours exclude larger couplings or mixing angles. Regions 
excluded due to perturbativity arguments are indicated; red, 
green and pink contours denote results from collider searches 
for invisible H decays [53], and mono jet [6] searches, and the 
LUX Collaboration m, respectively. 


assuming maximal coupling g x to dark matter. For the 
scalar-mediated model, limits are placed on the param¬ 
eter k x sin(# m i x ), where sin((9 m j x ) is the mixing angle 
between the scalar S boson and the Higgs boson, and n 
is a scaling constant; however, current calculations [2] 
of the gg —► HS production mode may be overestimated 
due to approximations made in evaluating the top-quark 
loop. 

In conclusion, a search for DM produced in association 
with a Higgs boson decaying to two photons has been 
conducted. Prior to these results, no bounds have been 
placed by collider experiments on the if+DM models 
discussed here. In addition, upper limits are placed on 
the cross section of events with large missing transverse 
momentum and a Higgs boson. 
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